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Abrikosov-Gor'kov (AG) theory, the foundation for understanding pair-breaking effects in con- 
ventional superconductors, is inadequate when there is an excitation gap (pseudogap) present at 
the onset of superconductivity. In this paper we present an extension of AG theory within two 
important, and diametrically opposite approaches to the cuprate pseudogap. The effects of impuri- 
ties on the pseudogap onset temperature T* and on T c , along with comparisons to experiment are 
addressed. 



Impurity effects in the high temperature superconduc- 
tors have been the subject of a large body of experi- 
mental and theoretical-literature concentrating on pair- 
breaking effects oiTgSJirHu'u, d-wave density of states ef- 
fects near .X, w 0,Btffl, local suppressions of the order 
parameter Jala, transport effectsj!3, r-and aspects of the 
superconductor-insulator transitiontil. Although there 
are works on the eflaaets of a single impurity in the_pfieu- 
dogap modelsE3E3oll3, with very few exceptionsQil3 lit- 
tle theoretical attention has been paid to the interplay 
between the widely observed cuprate pseudogap and the 
effects of disorder on pair-breaking. This is a particu- 
larly striking omission, given that j-a, major fraction of 
the superconducting phase diagramliZI is associated with 
a pseudogap. The goal of the present paper is to establish 
a formal mean-field structure (analogous to Abrikosov- 
Gor'kov (AG) theory) that incorporates this pseudogap 
in computing both T c and gap onset temperature, T*, 
along with other derived properties. Here we address 
two mean- field approaches ( orthogonal in their physics, 
but similar in their formalism), to the incorporation of 
the pseudogap: one in wiikhpthe pseudogap derives from 
superconductivity itselOcijEj ("intrinsic") and one in 
which it is "extpnsic" , either associated with a hidden 
order parametercM3, or .with band-structure effectsE3EJ. 
This intrinsic pseudogapuJo arises from a stronger than 
BCS attractive interaction which leads to finite momen- 
tum pair excitations of the normal state and condensate. 



In contrast to BCS theory, in the pseudogap phase 
there is an excitation gap present at T c , which, at low 
doping x, remains relatively T-independent for all T < 
T c e3. This necessarily will affect (i) fundamental charac- 
teristics of the superconducting phase as well as (ii) the 
nature of impurity pair- breaking. Indeed, to support (i), 
there are .strong indications from thermodynamics!!-] and 
tunnclingpil experiments that the-, effects of the normal 
state pseudogap persist below T<£A. Evidence in support 
of (ii) comes from the fact that pseudogap effects appear 
to correlate with the degree of the T c suppression in the 
presence of Zn impuritiesu. This suppression becomes 
progressively more rapid as the size of the pseudogap 
grows. 



Both intrinsic and extrinsic models for the pseudogap 
are associated with a generic set of mean-field equations. 
It is reasonable to stop at a mean-field level because 
these materials (in some, but not all respects) do not 
appear to be strikingly different from BCS superconduc- 
tors, and because the true critical regime appears to be 
rather narrowEI Moreover, we believe fluctuation effects 
around strict BCS theory such_as the phase fluctuation 
model of Emery and co-workerscJ are unlikely to explain 
the often very large separation observed between the gap 
onset temperature T* and T c . It seems more appropri- 
ate, thus to search for an improved mean field theoryEl 
Then additional fluctuation effects can be appended as 
needed. 

In this generalized mean field approach, in the clean 
limit and for T < T c , the gap and number equations are 
given by 
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where g sc is the coupling constant for the superconduct- 
ing order parameter, ip^ = (cos k x — cos k y ) is the e?-wave 
symmetry factor, A sc represents the superconducting or- 
der parameter, and A pg the pseudogap which persists in 
the T < T c phase. Finally, a is a ban d . in dex, which 
appears in some microscopic approachesEIlE^I to the ex- 
trinsic case. The momentum summation in the extrinsic 
case is over the reduced Brillouin zone. These equations 
depend in an important way on the electronic dispersion 
which differs in the two schemes. In the intrinsic school 
the fermionic dispersion is characterized by 
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Here £k is the "bare" band structure, taken to correspond 
to a nearest neighbor tight-binding model. This should 
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be contrasted with that in the extrinsic school, 
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The fermionic dispersions of the two schools differ as a 
direct consequence of the mechanisms that generate the 
respective pseudogaps. At the mean field level, a pseu- 
dogap due to pairing correlations forms as particles and 
holes mix to form the fermionic quasiparticles. Those 
of a spin- or charge- ordered state, though, are particle- 
particle mixtures. In the regime T < T c , where sharp 
excitations exist, these can be taken as the defining char- 
acteristics of "intrinsic" and "extrinsic" models of the 
pseudogap. Since generally /i ^ away from half-filling, 
only in the intrinsic school is the pseudogap pinned at 
the Fermi surface. 

The respective properties of the pseudogap lead to 
equations for its magnitude, which we summarize for T < 
T c , in terms of the particle-particle (x pp ) an d particle- 
hole (x ph ) susceptibilities. For the intrinsic school 
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Note that \ pp dependst3'EJ on the full excitation gap 
A. Here A pg (T) is associated with the number of finite 
momentum pair excitations of the condensate. These 
occur when the strength of the attractive interaction g sc 
is progressively increased, so that it is larger than that 
associated with the BCS regime. For the extrinsic school, 
the counterpart equation is 



X ^(0,0,A sc ,A ra ) = T 
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where g pg is the coupling constant for the pseudogap or- 
der and the momentum summation is over half of the 
Brillouin zone. Here we consider the pseudogap with 
same d-wave structure as the superconducting order. 

Figure [l] shows the temperature dependencies of the 
different energy gaps obtained by solving the complete 
set of equations in the two pseudogap schools within the 
underdoped regime. In the intrinsic case T* marks the 
gradual onset of the pseudogap, which is associated with 
bosonic or pair excitations formed in the presence of a 
stronger-than-BCS attractive interaction. Only at and 
below T c does the identification of A become precise, so 
that for this (intrinsic) case we plot an extrapolation of 
Eqs. 00 and | to T > T c . Figure 0a shows that be- 
low T c the fraction of the bosonic population joining the 
condensate of zero-momentum pairs (oc A^ c ) increases 
at the expense of the finite-momentum bosonic fraction 
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FIG. 1: Energy gaps for intrinsic (a) and extrinsic (b) cases. 
Solid lines are the total excitation gaps, dotted lines the super- 
conducting order parameters and dashed lines the pseudogaps 
below T c . The gaps are in units of 4t^ . The curve for T > T c 
in (a) represents a rough extrapolation. 



(oc Ap g ) until the fully condensed ground state is reached. 
By contrast, for the extrinsic case (Fig. |l]b) superconduc- 
tivity forms on top of a pre-existing excitation gap in the 
effective band structure which first appears at T*, the 
phase transition temperature marking the onset of the 
extrinsic order. 

One can capture the key physics of these two schemes 
in a reasonably accurate phenomenological approach. 
The bpaonic excitations associated with the mean-field 
theoryL3 of Eqs. (la), ([it]) and (||) lead to the tempera- 
ture dependence of the pseudogap below the clean limit 
critical temperature T c q 
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These bosons are, thus, associated with a quasi-ideal 
Bose gas. By contrast for the extrinsic case, in the well- 
established pseudogap regime, below T c q, the pseudogap 
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FIG. 2: DOS for intrinsic and extrinsic models at T = 0. 
Only one gap structure appears in the intrinsic DOS, while 
two distinct gap structures appear in the extrinsic DOS. 



is relatively T-independent 



A 2 p JT) = A 2 (T c0 ), T<T CI 
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Here we define A(T c0 ) = A pg (T c0 ). In both Eqs. (§) and 
(0) above, we may view A(T c o) as, a phenomenological 
parameter taken from experiment!!-! We will adopt this 
approach here in large part because it provides a more 
readily accessible theoretical framework for the commu- 
nity, and because it connects more directly with experi- 
ment. 

The pronounced differences between the fcrmionic dis- 
persion in these two theoretical schools can be seen from 
the associated densities of states (DOS) plotted in Fig. [|, 
which compares the intrinsic and extrinsic models at 
T = 0. In the intrinsic model one sees only one exci- 
tation gap featureE] 
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in Fig. |[ cen- 
tered around the Fermi energy. Van Hove singularities 
are also apparent here as relatively sharp structures. In 
contrast, there exist two distinct features for the extrin- 
sic theory. The more prominent pseudogap peaks are 
centered around — /i, while the .superconducting peaks 
appear around the Fermi energyLj. Indeed, for this ex- 
trinsic case, only in the limit /i = can one readily define 
an excitation gap A as in a conventional superconducting 
phased, satisfying Eq. (2b). That the superconducting 
order parameter and pseudogap contribute to separate 
features in the density of states represents a rather clear 
signature of this extrinsic pseudogap school. To date, the 
bulk of experimental tunneling data supports a-»ieturc 
in which there is a single excitation gap feature! 3 ni 33 !, al- 
though there are some reports of multiple-gap structures 
in c-axis intrinsic tunneling spectroscopy 3 .! At T = T c , 
the extrinsic superconducting gap closes and the densi- 
ties of states for the two schools become quite similar, 
save for the pinning of the gap minimum to the Fermi 
surface in the intrinsic case. 

We turn now to impurity effects which, just as in the 
BCS case, are not expected to change the formal struc- 



ture of our mean field theory. The greatest complica- 
tion is associated with the impurity- renormalized A pg , 
calculated from all possible diagrammatic insertions of 
the impurity vertex into the particle-hole and particle- 
particle susceptibilities [see Eqs. (|4|) and (JsJ) .] A detailed 
study of these effects in the intrinsic case appears in 
Ref. [l6| although here we will proceed more phenomeno- 
logically within both schools. We base the present treat- 
ment on analogs of the clean limit mean field gap equa- 
tions Eqs. (JTa|) and (lb) with substitutions A sc — > A sc , 



A pg — > A pg , oj n — > u> n , and fj, — > p,. At the phenomeno- 
logical level the T-dependence of the intrinsic pseudogap 
is given by 



A 2 (T)«A 2 (T C ) 



and for the extrinsic case, 
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where, in both schools, the excitation gap A(T C ), is pre- 
sumed to be determined from experiment. 

To compute the renormalized frequency iLu(iuj n ) and 
chemical potential fi(iu) n ), we follow the usual im- 
purity T-matrix approach. We presume an s-wave 
short-range impurity potential V(r) = uS(r — r,). 
The impurity scattering matrix T(uj n ) in Nambu space 
satisfies the Lippman-Schwinger equation: T(uj n ) = 
u&3 ^1 + T(uj n ) <?(k, uj n )^j , where g is the impurity- 
dressed Green's function, 
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Here A is either the full gap or superconducting order 
parameter in the intrinsic and extrinsic cases, respec- 
tively, and &i are Pauli matrices. Here we suppress the 
band index in the extrinsic case. Labeling components 
as 9 = Si 9i®ii the regular and anomalous Green's func- 
tions are G = go + .93, F = —g\. The frequency and 
chemical potential are renormalized through impurity 
self-energy E = n;T, and ioj n — iui n — So, p, = fx — S3, 
where rii is the number of impurities per unit cell. We 
note that the T-matrix for the extrinsic school depends 
only on the band structure and is independent of the 
specific type of extrinsic order. 

The components of the self-energy are given by 
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There is no frequency-dependent self-energy associated 
with gap renormalization due to d-wave symmetry. Fi- 
nally, the magnitudes of A, A sc and A pg can be obtained 



FIG. 3: Intrinsic DOS for the clean and dirty cases at T = 0. 
The DOS is centered around the Fermi energy. 
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FIG. 4: Extrinsic DOS for the clean and dirty cases at T — 0. 
The DOS is centered around — fi. 



using Eqs. (@),(||) and (||), presuming that the excitation 
gap at T c is taken from experiment. Here we take the bare 
lattice dispersion £k = — 2t» (cos k x + cos k y ) — 2t± cos k z 
so that the dimensionless coupling constant is given by 
<?/4*||. 

Figures |3| and || show the effects of impurities (for 
unitary scattering) on the density of states at T = 0, 
in the intrinsic and extrinsic cases respectively. As can 
be seen, particularly for the intrinsic case, impurities de- 
crease slightly the height and separation of the gap peaks 
[See Fig. |5j below] and fill in the low frequency region, but 
otherwise their effects are not dramatic. For the extrin- 
sic school, the superconducting gap region is more qual- 
itatively affected by pair-breaking, while the pseudogap 
peaks remain relatively robust. It can be inferred from 
these figures that with increasing disorder the differences 
in the two schools diminish, from the perspective of the 
density of states, except that the position of the mini- 
mum in the extrinsic case is not tied to the Fermi energy. 
Physical differences, however, remain profound, particu- 



FIG. 5: Temperature dependencies of the full excitation gaps 
for the intrinsic case at strong (<?/4tii = — 1.2) and weak cou- 
pling (<?/4t|i = —0.15, inset), in the unitary limit at different 
impurity densities. Temperatures are normalized to the clean 
limit Tq and gaps are normalized to the zero-temperature 
values in the clean limit Aqo- 



larly in the electrodynamicsL^ES of the superconducting 
phase. 

In the remainder of this paper we focus on the behavior 
of T c (and T*) and the appropriate generalization of AG 
theory in the presence of a pseudogap. For definitcness, 
we consider Zn doping experiments where we exploit the 
experimental observation that the excitation gap A at 
T c is relatively unaffected by Zn impurities^. We focus 
here on the unitary limit (1/u = 0), which is regarded as 
relevant to Zn doping in the cuprates. 

We begin with the intrinsic school, where the sensitiv- 
ity of T c and T* to the impurity concentration rij can 
be studied as a function of a single coupling constant 
g = g sc , which we presume to be unaffected by the addi- 
tion of impurities. Figure]^ shows the behavior of the ex- 
citation gap A(T) vs temperature normalized to its clean 
limit value, obtained using the impurity-generalized form 
of Eqs. ( p"a| ) and (lb). The figure should be viewed as 
extending above T c only-in the sense that it provides a 
reasonable extrapolation^!] as well as estimate of T* . In 
reality, Fig. 0a indicates that a crossover description for 
the excitation gap at T above T c is more correct. The 
main panel corresponds to the strong (g/At\\ = —1-2) and 
the inset the weak (g/At^ = —0.15) coupling regimes for 
various values of the impurity density rij in the unitary 
limit. In the weak coupling regime, the suppression of the 
gap is largest, as is the reduction in T*. In the strong 
coupling case, the suppression is smaller and at low impu- 
rity densities, the curves are very close to those obtained 
in the clean limit, indicating smaller pair-breaking effects 
on the excitation gap and its onset temperature T*. 

Figure ^ shows the way in which impurities suppress 
the phase coherence temperature T c at different coupling 
strengths (in the unitary limit), based on the assump- 
tion, supported experimentallyu, that the excitation gap 
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FIG. 6: T c suppression due to impurities for l/u = (unitary 
limit) in the intrinsic case. The temperatures are normalized 
to the clean limit T c o- 
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FIG. 7: Initial slopes of T* and T c suppression vs 
doping, in the unitary limit for the intrinsic case. The in- 
set presents a comparison between theory and experimental 
data from Ref. |§| See text for details. 



at the appropriate T c is relatively independent of impu- 
rity concentration. It can be seen that the suppression 
rate increases as the coupling becomes stronger, or effec- 
tively as A(T C ) increases. Similar results for the extrinsic 
case were obtained in Ref. ||. This faster T c suppression 
in the strong coupling regime can be understood through 
a simple physical picture. Impurity scattering will pro- 
duce states which fill in the gap and eventually destroy 
superconducting coherence. In the strong coupling (pseu- 
dogap) regime, where the normal state already has a gap, 
fewer impurities are required to restore the system to the 
"normal" state. 

We turn now to calculations which can be directly com- 
pared with experiment and plot the normalized slopes of 
T* and T c with respect to increasing Zn concentration, 
for varying hole concentration x, first for the intrinsic 
case. To convert from the coupling constant parameter g 
to x we take as input the experimentally measured val- 
ues of p s (x,0) and the measured excitation gap at T c . 
Here it is adequate to choose these values corresponding 
to the pristine case, and presume that Zn doping does 
not affect the excitation gap at T c . Figure indicates 
the initial slope ( \/T$dT /drii, where Tq is the appropri- 
ate clean limit temperature) for T* (dashed line) or T c 
(solid line). In the overdoped limit, the theory is asymp- 
totically equivalent to standard AG theory, in which also 
T* = T c . For smaller values of x the slope decreases so 
that T* is only weakly dependent on impurity concentra- 
tion. By contrast, the initial T c slope (solid line) shows a 
very different hole concentration dependence. As the hole 
concentration decreases, the slope decreases. However, 
in the very underdoped regime, where the pseudogap is 
well established, the curve turns around and rapidly in- 
creases. The inset presents a comparison of theory and 
experiment^ as r\ — (dT c /dni)/(dT c /drii,x — 0.20) vs 
z = A pg (T c )/(Apg(T c ),x — 0.05), where the agreement 
appears to be reasonable. There are fewer systematic 
studies of impurity- induced changes in T*; however, the 
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Initial slopes of T* and T c suppression -^r vs dop- 



ing, in the unitary limit for the extrinsic case. The inset 
shows the clean phase digram given by coupling constants 
g pg /4t|| = —0.4 and g ac /At\\ = —0.375. The critical doping x c 



where A P9 vanishes is around 0.15. 



small effect found here at low x appears to be compatible 
with the data. 

Finally in Fig. || we present the counterpart plots of 
the initial slopes for T c and T* in the extrinsic case. T c 
is computed as in the intrinsic case by assuming A pg (T c ) 
is relatively insensitive to impurities. Impurity renor- 
malizations are determined through Eqs. ( pd| ) and (12) 
while the suppression of T* is calculated via Eqs. (lb) 
and (H) , extended to include appropriate impurity renor- 
malizations. The inset shows the clean phase diagram 
which forms the basis-,for these calculations. Our fit 
to the published forn£3 of this phase diagram provided 
values for the coupling constants gpg/^t\\ = —0.4 and 
g sc /At^ = —0.375. To make contact with experiment 
we chose a parameter set in which T c /T* and n s /A(0) 
were reasonably well fit to experiment in the underdoped 
regime. As is similar to the intrinsic case, there is a 
dramatic increase in the slope of T c as the insulator is 
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approached. This increase is associated with the onset 
of the pseudogap which occurs for x < 0.15. Above this 
critical concentration T* is zero, and the system becomes 
a conventional dirty BCS superconductor. In this way, 
the intrinsic and extrinsic schools differ, since for the for- 
mer at large x, T* — > T c . 

The theoretical machinery that we have set up has 
strong similarities-, to an approach taken by Loram 
andjjjollaboratorsa, extended further to the disordered 
caseoQ. It should be stressed, though, that their ap- 
proach is a hybrid of extrinsic and intrinsic pseudogap 
theories, where the temperature dependence of the var- 
ious gap parameters corresponds to the extrinsic case 
(shown in Fig. |l|b), whereas the dispersion and superfluid 
density corresponds to an intrinsic pseudogap. As shown 
in this paper, pair-breaking effects on T c can be success- 
fully addressed at a semi-apantitative level both in in- 
trinsic and extrinsic modelsu. It should be noted that the 
rather strikingly different sensitivities of T* and T c to im- 
purity concentration which are found experimentally, are 
often taken as an indication that the cuprate pseudogap 
cannot be intrinsic, i.e., related to the superconductivity, 
itself. Indeed similar results are found in the presence of 
magnetic field pair-breakingES and it should be viewed 
as one of the fundamental results of this paper that this 
inference is incorrect. The differences lie in the fact that 
an excitation gap is present when T c is established, but 
not at T* , and it is this gap in the density of states 
that contributes to the stronger pair-breaking effects in 
T c . Indeed, it is precisely this excitation gap which in- 
validates the results of conventional AG theory. It may 
be necessary eventually to incorporate an even more local 
treatment of pair-breaking than that discussed here, but 
such a Bogoliubov-de Gennes generalization must include 
pseudogap effects. Indeed, the very basis for a more local 
treatment of impuritiesa is the observed small coherence 
lengths, which are at the heart of the present "intrinsic" 
pseudogap theories^. 



In summary, in this paper we find within two dia- 
metrically opposed pseudogap schools, that pseudogap 
effects at and below T c must play an essential role in 
pair-breaking. While there is no definitive experiment 



to distinguish 
gued elsewhereE 



reen these two schools, we have ar- 
3 that the intrinsic dispersion leads 
to smaller and more benign modifications of BCS the- 
ory. In both theoretical approaches, the rather robust 
behavior for T* and the associated excitation gap in 
the underdoped regime, found in the presence of im- 
purities may be associated with the widely observed 
superconductor-insulator transitiontil. Superconducting 
coherence is more readily destroyed than is the excita- 
tion gap (and T*), thereby leading to an insulating state 
when T c is suppressed to zero, in much the same way as 
in the presence of applied magnetic fieldsc3. While there 
are clear differences, seen particularly in electrodynami- 
cal calculationst§E3 (as well as density of states effects) 
between the intrinsic and extrinsic pseudogap schools, 
the impurity sensitivities of T c within these two different 
approaches are quite similar, and reasonably consistent 
with experiment. This similarity derives from the fact 
that both mean field theoretic calculations of T c have 
a general BCS-like character, except for the presence of 
a (pseudo) gap at the onset of superconductivity. For 
T* the differences are more apparent in the overdoped 
regime and this, in turn, reflects the fact that T* — » in 
one case (extrinsic) , whereas T* — > T c in another (intrin- 
sic). In this paper we have set the stage for a computa- 
tion of transport properties which require as an essential 
input, an understanding of impurity effects. The general- 
ization of AG theory presented here should help to clarify 
the important role played by pseudogap effects, at T c and 
their relation to impurity-induced pair-breaking. 
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